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Abstract 

Measurements of the B° S -B Q S mixing phase <f> s appear to be converging towards the SM prediction, suggesting that 
contributions from New Physics (NP), if present, are small. This poses the question of whether smallish NP in <p s can 
be distinguished from the hadronic uncertainties present in the interfering B s decay mode. In this paper we discuss 
the potential of extracting <p s from the decay modes B s — > Jjif/rf'^ and B s — > 7/(^/0(980) and how their respective 
hadronic uncertainties can be controlled. In addition, we demonstrate how the branching ratios of the former decays 
can be used to estimate the 77-77' mixing angles. Finally, we point out that effective lifetime measurements of decay 
modes such as the B s — > //i/r/o(980) can be used to constrain both <f> s and the B s decay width difference Ar s , 
complementary to the usual time-dependent tagged analysis. 
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1. Introduction 

A current trend in flavour physics is that CP observ- 
ables are approaching their Standard Model (SM) pre- 
dictions. A notable example is the B^—B® mixing phase 
(f> s , which is expected to be very sensitive to new CP vi- 
olating contributions (for a review see Ref. HI). Specif- 
ically, in the SM this phase is predicted to take the tiny 
value <^ M = -2.1°, towards which the latest experimen- 
tal results from the flagship B s — > J/i//<p analysis now 
appear to be converging J31 : 



(0.0 + 6.0)° 

(-31-5!^)° 
[-34.4,6.9]° 



LHCb 0, 

D0 0, (1) 

CDF 0. 



Included for correctness is a possible hadronic phase 
shift A0j^ . . for each of the final state polarizations A (6). 
As these shifts are expected to be small, we conclude 
that there is no large smoking gun signal of New Physics 
(NP) in B°-B° mixing. 
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However, as precision improves, what will we be able 
to conclude about smallish NP in <p s l To distinguish 
it from SM physics, it will become necessary to con- 
strain the theoretical uncertainties in parameters such as 
A</>j^ to a sufficient accuracy. Because these uncer- 
tainties are driven by the hadronic nature of these de- 
cays, this is not an easy task. Standard tools include 
lattice QCD, QCD sum rules and SU(3)p flavour sym- 
metry relations between decays. Concurrently, it would 
be helpful to overconstrain the phase <p s with different 
analyses, analagously to the determination of the angles 
of the Unitarity Triangle. In this paper we explore de- 
cay modes and analysis strategies complementary to the 
usual time-dependent, tagged B s — > J/if/</> analysis. 

The layout is as follows. We first consider the op- 
timal decay mode structure B s — > J/iffss and how its 
hadronic uncertainties can be controlled. We then dis- 
cuss, in turn, two alternative final states with this struc- 
ture: B s -> J/if/rf^ and B s -> 7/^/ (980). Finally, we 
highlight the utility of effective lifetime measurements, 
observables that do not rely on flavour tagging, for con- 
straining the phase <f> s and the B s decay width difference 
AF S . 
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2. CP violation in B s -* J /if/ ss decay modes 



The standard manner of measuring the CP violat- 
ing phase <p s present in the B° s -B® mixing amplitude 
is via its interference with a specific decay amplitude 
for a final state /. In this paper we consider only CP- 
eigenstates and denote the CP eigenvalue by r\f. Be- 
cause we measure the interference between mixing and 
decay, the CP violation of the decay mode B s — > / is 
also relevant and can be parameterised by the direct CP 
asymmetry Cf and a phase shift A<f>f. In practice what is 
measured is the time-dependent rate asymmetry, which 
depends on the parameters <f> s , Cf and A(f>f as follows: 

r(B°(o^/)-r(5?(?)^/) 



B, 



r(ii?(f)->/)- 1X32(0 -»/) 

C f cos(AM s t)-S f sin(AM. s f) 
cosh(Ar s f) + J{{ r sinh(Ar j t) 



(2) 



with 



M f &r + iS f = -rj^\ -C}exp [i(<f> s + A0y)] . (3) 

The interfering phase is therefore the combination <p s + 
A(f>f and the extraction of <f> s is limited by any theoretical 
uncertainty on the phase A(pf. 

An optimal choice for the decay mode of a B s inter- 
ference measurement is a B s — > J/if/ss transition, where 
the ss quark pair forms a bound state. This is because all 
the topologies contributing to this decay that contain CP 
violating phases are doubly Cabibbo suppressed. Let- 
ting e = A 2 /(I - A 2 ) ^5% parameterise the Cabibbo 
suppression and explicitly extracting the CP violating 
phases, the B s — > J/if/ss decay amplitude in the SM has 
the structure: 



M(B S -> J /if/ ss) oc 1 + ee' y b e" 



(4) 



where y is the usual angle of the CKM Unitarity Tri- 
angle. The complex parameter b e' 9 is the ratio of CP 
violating over CP conserving decay topologies. For the 
transition in question the CP violating topologies are all 
penguin in nature whereas the CP conserving topolo- 
gies include, along with penguins, a colour-suppressed 
tree topology. As the individual decay topologies de- 
scribe hadronic, non-perturbative, processes, the param- 
eters b and 9 are difficult to compute theoretically. For- 
tunately, assuming the tree topology is dominant, the e 
suppression ensures that the CP violation present in the 



J/if/ss is small at best: 

C/M S s = -2eb sin6» sin y + 0(e 2 } 



-5% x 



siny 



sin 68° 



0.5 



sin 9, 



j/^ss -2eb cos 9 sin y + O (e 2 ) 



3°x 



siny 



sin 68° 



0.5 



cos 9. 



(5) 



(6) 



However, with the increasing level of precision to be of- 
fered by the LHC, these effects may no longer be negli- 
gible and it will be important to constrain the underlying 
hadronic physics. 

A standard tool to control the hadronic parameters b 
and 9 is via the SU(3)f flavour symmetry of strong in- 
teractions. In the limit that the three lightest quarks have 
masses that are degenerate and much lower than the 
Aqcd scale, the strong dynamics of the various topolo- 
gies is invariant under their interchange. Consequently, 
in this limit, the hadronic parameters b and 9 are also 
common to B q — > J/ifi dq decays for q e {u,d, s}. The 
key difference of the latter decays is that the CP violat- 
ing decay topologies are not doubly Cabibbo suppressed 
i.e. the parameter b is no longer accompanied by a factor 
of e: 

M(B q -> J /if/ dq) SU £ 3)f 1 - e iy b e ie , (7) 



These decays thereby act as control channels, allow- 
ing us to constrain the parameters b and 9 up to SU(3)f 
breaking corrections. 

A prime example is the golden decay mode Bj — > 
J/i//Ks, which has been controlled using the decay 
B c i — > J/iffx° to constrain the hadronic phase shift 
A<f to an accuracy of ±1.5° 0- In the fu- 
ture the decay mode B s — > J/if/K s should improve on 

this num. 

The question now is what to choose for the ss me- 
son in the B s — > J /if/ ss transition. The vector meson (f> 
is believed to be almost completely ss and is what was 
chosen for the flagship analysis B s — > J/if/<p. However, 
this results in a vector-vector final state that requires a 
complicated time-dependent angular analysis in order to 
resolve the CP violating phase </> s + M^jm introduced 
in ([T]i ifTTl fl2ll . In Ref. |6j the hadronic control chan- 
nels are studied for the B s — > J/tf/(f> decay. Other bound 
states that are believed to contain an ss component are 
the pseudo-scalars 77 and 77' and the scalar /o(980). Their 
key advantage is that they do not require an angular 
analysis. In the sections that follow we discuss their 
feasibility as alternative final states. For a general dis- 
cussion of B s — > J /if/ ss modes see also Ref. 1 13 1. 
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3. B s -> Jlilftf* 

The pseudo-scalars rj and 77' are believed to be a mix- 
ture of isospin singlet states: 

\q) = cos (pp (\uu) + \dd)\ - sin^p|ss) 

V2 v ' 

\rf) = cos <pc sin (pp — (\uu) + \dd)) 
V2 v ; 



• cos (p c cos (p P \ss) + sin (f> G \gg), 



(8) 



The mixing angles are not yet tightly constrained, with 
30° < (/>p < 45° and the possible glueball state given 
by \<pc\ ~ 20° (for details and a review see Ref. [14|). 
Because the final state is no longer purely J/ipss but 
now also contains J/ip{uu,dd) contributions dependent 
on the strength of the mixing angles, there are additional 
exchange and penguin annihilation topologies present. 
Fortunately, in the SM all CP violating topologies are 
still doubly Cabibbo suppressed for the final state. 

TheB s — > J/xprf'^ decay modes are related by SU(3)p 
flavour symmetry to the B d — > I/if/rf^ modes, respec- 
tively. In Ref. |fl5i this symmetry relation was em- 
ployed in a proposed strategy to determine the CKM 
angle y. Instead, assuming we have measured y, the 
Bd — * •> 70 "H W decay modes can act as control chan- 
nels for an extraction of the mixing phase <p s from 
B s — > J/tf/rf') decays, as discussed in Sectional 16|. 

Due to the presence of the mixing angles, however, 
there is in general no clean way to relate the hadronic 
parameters b and 9 of the former decays to the latter in 
the limit of SU(3)p flavour symmetry. A relation can in 
general be made if we may assume a hierarchy between 
topologies. Namely, if we may assume that the strong 
amplitudes of the tree and penguin topologies contribute 
much more than the exchange and penguin annihilation 
amplitudes. Such hierarchies can be determined by as- 
suming flavour symmetry and considering bounds on 
branching ratios of decays where exchange or annihi- 
lation amplitudes are dominant. For example, a bound 
on the CP violating exchange topology could be derived 
from an updated B s — > J/if/n branching fraction, the 
current value of which dates back to an L3 measured 
from 1997 IT71. 

In terms of experimental progress, the Belle collabo- 
ration has performed the most accurate branching ratio 
measurements lfl8l[T9ll : 

BR(B S -> J/t/nf) = 5.10!,^° x 1(T 4 , 
BR(B S J/i//T]') = 3.71 + ';°° x 1(T 4 , 

BR(B d J/iJ/tj) = 12.3+J;! x 1(r& > 
BR(B d -> J/iffT]') < 7.4 x 1(T 6 . 
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Figure 1 : Constraints on the rj-rf mixing parameters from the B s ^ — » 
and B d — > J/ij/n branching ratios, as discussed in the text. 
Note that not all discrete amibiguities following from the constraints 
in j9) are shown. 



However, the prospects of measuring the CP observ- 
ables with a time-dependent tagged analysis at the LHC 
are not good. This is because 77 and 77' decay promi- 
nently to photons and neutral pions, which are difficult 
experimental signatures for the LHC. Measurements be- 
yond branching ratios will likely have to wait for the 
Super "flavour" factories. 

If we assume the topological hierarchy from above, 
it is possible to use the branching ratio measurements 
to constrain the 77-77' mixing angles. To this end, we 
construct the ratios: 



R, = 



BR{B S -> Jjftf) ( <D 
BR(B S 

COS 2 (fir 



R = 



tan z cpp 
BR(B d 



0.91+0.18 



Jim 



BR(B d -> J/tpn°) 
cos 2 P = 0.77 ±0.14 



4>", 

a 



BR(B d -> JIM) 



BR(B d 



J/lffTj) 



tan" i 



(9) 



where is the phase space factor for a/3, — > J/yj P 
transition. Similar expressions ignoring a possible glue- 
ball component ((pc = 0) are given in Ref. EH . In Fig- 
ure [T] we show the corresponding curves on the (p<j-<pp 
plane. The contour for R d is included for illustration and 
has been engineered to meet at the intersection point of 
the other two contours. Due to its sharp intersection 
with the R s contour, future measurements should allow 
it to strongly constrain the mixing angles. 
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4. b s jwum) 

The decay B s — > J/iftfi)(980) was originally proposed 
as an alternative to B s — > in Refs E2l l23l . Al- 

though its branching fraction, with /o(980) — > n + n~, 
is approximately a fourth that of B s — > J/i[t(p, with 
— > K + K~, its chief advantage is that it does not re- 
quire an angular analysis. Due to the charged particles 
in the final state, this decay mode, unlike B s — > Jlijjrf'\ 
is experimentally viable at the LHC detectors. In fact, 
the LHCb experiment has reported the following mea- 
surement of its CP violating phase [24 1 

^ + A0 Wo = (-25±25)°. (10) 

Strictly speaking, the experiment has ignored the pos- 
sible hadronic phase shift in the decay mode. As we 
will now argue, the theoretical uncertainty on A<pj/^,f is 
significant and should be included in the future |251 . 

The underlying source of uncertainty in the CP vio- 
lation present in the B s — > //i^/o(980) transition is the 
uncertain nature of the /o(980). There is no consensus 
in the literature on what exactly this scalar meson is and 
how much of it is dominated by a ss combination. Two 
popular pictures of the /o(980) are that it is a conven- 
tional qq scalar meson or a tetraquark. More exotic in- 
terpretations include a KK molecule (for a review see 
Ref. ||26| and references therein). 

In the qq description, the /o(980), together with the 
cr(600), can be a mixture of isospin-singlet states with a 
mixing angle ip M in analogy to the angle <p P for the rf'^ 
states defined in (|S). The puzzling feature of this pic- 
ture, however, is why the scalars, with nonzero angular 
momentum, are lighter than their pseudo-scalar coun- 
terparts. 

In the tetraquark picture, the /q(980) is a bound state 
of a diquark and an anti-diquark in a low energy con- 
figuration with no orbital angular momentum. This in- 
terpretation reproduces the mass ordering of the light 
scalars in the expected way, and a mixing of less than 5° 
is expected with the cr(600) tetraquark state (which we 
ignore in our analysis). In Ref. 11271 a consistent picture 
of all the scalar mesons is presented, where tetraquark 
and qq states mix via instanton effects. A fit to the mass 
spectrum results in the light scalar mesons, with the 
/o(980) among them, being predominantly tetraquark. 

Which topologies contribute to the hadronic parame- 
ters b and 9 is strongly dependent on the /o(980) picture 
that one subscribes to. In the conventional qq picture the 
distribution of topologies is related to the mixing angle 
(pM similarly to the J/tfrq^ final state. In the tetraquark 
picture an entirely new topology is possible, which we 



4 



b 




Figure 2: The additional topology A± q that contributes to the B s — > 
J/ifr fo(?SO) transition in the tetraquark picture of the /b(980). 

label A4 q and depict in Figure [2] whose relative hierar- 
chy to the other topologies is not known. 

A potential control channel for constraining the 
hadronic parameters b and 9 is the transition Bj — > 
/o(980). However, in both pictures the dependence 
of the respective parameters on the various topologies 
in the SU(3)p limit is sufficiently complex that we must 
assume a hierarchy for the topologies. Namely, that the 
tree and penguin amplitudes dominate over the rest. A 
possible channel to constrain the new topology A\ q is 
B s -> J/i//a°(9S0) in the tetraquark picture (25). 

To estimate the current uncertainty on the CP viola- 
tion in the B s — > 7/i/'/o(980) decay mode, we make 
the conservative assumption that the tree amplitude is 
greater than or equal to the other topologies and leave 
the strong phase 9 unconstrained. Setting y = (68 ± 
7)° ED, it then follows from ^ that 

A0 Wo € [-3°,3°], \C mfo \ < 0.05. (11) 

Combining this with the SM prediction for the CP vio- 
lating phase and inserting it into (|3]l gives the SM 
prediction for the mixing-induced CP asymmetry: 

5 Wo | sm 6 [-0.09,-0.01]. (12) 

This should be compared to the naive SM estimate 
sin^IsM = -0.036 ± 0.002. 

The SM prediction for the branching fraction of our 
proposed control channel is 

BR(B d ^J/if,f {) ;fo->7T+n-) 

If tatupM 1^ 
Ltan35°J • 99 ( 13 ) 

1 : tetraquark 

where a mixing angle of tpM = 35° in the conventional 
qq picture gives an equivalent rate to the tetraquark pic- 
ture. Thus an observation of this decay would signal 
either a significant dd component for the /o(980) or hint 
at its tetraquark nature. 
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5. B°-B° mixing constraints from effective lifetimes 

Instead of performing a full fit to the tagged time- 
dependent rate for a B s — > / transition as given in Q, it 
is also possible to construct observables sensitive to CP 
violation from the time-dependent untagged rate: 

<r» = r(B,(o -» f) + r(B s (t) -» /) 

oc <r' /T «< [cosh( Ar, f) + yi{ r sinh(AF. s t )] , (14) 

where Ar s = F^ - Fj^ is the B s decay width difference 
and tb % the mean B s lifetime. Given enough statistics, 
it is possible to fit for both AF, and J^ hT separately. 
Alternatively, given limited statistics, one can instead 
fit a single exponential, which yields the effective life- 
time. Fhe effective lifetime obtained from a maximum 
likelihood fit to the untagged rate is analytically equiv- 
alent to the time expectation value (see the appendix of 
Ref. [29 1). We may therefore define the effective life- 
time as oni 



Teff = 



r t(r f )dt 



Tb 



1 



■ v v 



(15) 



where y s = t^AF.s/2. Fhis observable gives us access 

-A<P f 
also 



to yv kY defined in (F3J, which in turn depends on < 
for the given final state /. This dependence on Ji J hT , 
makes it a key ingredient for correcting experimentally 
measured branching ratios in the presence of a nonzero 
decay with difference [29 1. 

The most accurate determination of the B s — > 
//^/ (980) effective lifetime is by LHCb ED: 



r Wfo = 1.70 + 0.040 ± 0.026 ps. 



(16) 



Interestingly, this measurement deviates by more than 
1 <x from the SM expectation: 



r Wo sM = (1.616 ± 0.034) ps, 



(17) 



obtained from the theoretical estimate ArJ h = 0.087 ± 
0.021 ps" 1 021 and th = 1.509 ± 0.017 ps 03). 

To avoid using theoretical input for AF S , we can in- 
stead treat it, along with <f> s , as one of the parameters we 
seek to determine. As we have theoretical estimates for 
Cj/^/ , A</>j/j,f and a measurement for T//^/ , we are left 
with just these two unknowns. The expression in ( fT5j ), 
which is a cubic equation in AT S , can consequently give 
a contour in the <f> s - AT S plane. Although this by itself 




TJIth = [1-700 ± 0.040 ± 0.026] ] 
1 tk+k- = [1.455 ± 0.046 ± 0.006] ] 
Z2 68% CL of x 2 & 



-90 -45 45 90 

<k Peg] 



435 



Figure 3: The measurements of the effective B, — > J/ifr /o(980) and 
B s — > K + K~ lifetimes projected onto the (p s -AT s plane. The shaded 
bands give the 1 cr uncertainties of the effective lifetimes. The 68% 
confidence region originates from &X 1 fit 1281 . 



is already useful, the expression in ( fT5] l has another in- 
teresting feature. Consider the decays of a B s meson 
to a pair of CP eigenstate final states, one CP even and 
one CP odd, with a comparable CP violating phase in 
the transition. Due to the dependence of Jl Ar on rjf, as 
given in Q, their respective effective lifetimes will give 
distinctly behaving contours in the <f> s - AT S plane. The 
intersection of these contours then pinpoints <f> s and AT S . 

Auspiciously, the only other effective lifetime cur- 
rently measured is for the B s — > K + K~ decay mode, a 
CP even final state expected to have a small CP violating 
phase. LHCb has the most accurate determination (34): 



Tk+k - = 1.455 ± 0.046 ± 0.006 ps. 



(18) 



In a separate flavour symmetry analysis the CP violation 
of this decay was estimated using B factory data from 
its t/-spin (a subgroup of SU(3) F ) partner decay B d — > 



n n , giving 



%r =-(l0.5!|9°, C Pr =0.09. 



(19) 



In Figure[3]we show the two contours for the CP even 
and odd final states on the (f> s -AT s plane. The 1 cr error 
bands correspond to the errors on the effective lifetime 
and the 68% confidence region is the result of a com- 
bined x 1 fit- The individual 1 cr confidence level results 
of the fit for the solution closest to the SM are 



<p s = (-49!J?)° , AF. S = 0.189!°;°g ps - 



(20) 



which, interestingly, do not agree well with <p s = (0.0 + 
6.0)° and Ar. s = (0.116 ± 0.019) ps 1 coming from the 
flavour tagged time-dependent LHCb B s — > J/ifr(p anal- 
ysis Q. The <p s value does agree with the less precise 
D0 and CDF results given in ([TJ. 
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A clear advantage of the effective lifetime strategy 
is that it reveals the behaviour of the untagged data. 
This information can be hidden in the timed-dependent 
tagged analysis given in Q, where the tagged and un- 
tagged data are generally combined into a single fit. 
Having multiple distinct contours pinpointing a solu- 
tion in the (f> s ~£T s plane would complement informa- 
tion from a global fit in a useful way, in analogy to the 
famous Unitarity Triangle plot. 
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